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Abstract 

This work is devoted to the study of the decay of multiscale deterministic solutions 
of the unforced Burgers' equation in the limit of vanishing viscosity. 

It is well known that Burgers turbulence with a power law energy spectrum Eo(k) ~ 
\k\ n has a self-similar regime of evolution. For n < 1 this regime is characterised by 
an integral scale L(t) ~ t 2 /( 3 + n ) ) which increases with the time due to the multiple 
mergings of the shocks, and therefore, the energy of a random wave decays more slowly 
than the energy of a periodic signal. 

In this paper a deterministic model of turbulence-like evolution is considered. We 
construct the initial perturbation as a piecewise linear analog of the Weierstrass func- 
tion. The wavenumbers of this function form a "Weierstrass spectrum", which ac- 
cumulates at the origin in geometric progression. "Reverse" sawtooth functions with 
negative initial slope are used in this series as basic functions, while their amplitudes 
are chosen by the condition that the distribution of energy over exponential intervals 
of wavenumbers is the same as for the continuous spectrum in Burgers turbulence. 
Combining these two ideas allows us to obtain an exact analytical solution for the 



velocity field. We also notice that such multiscale waves may be constructed for mul- 
tidimensional Burgers' equation. 

This solution has scaling exponent h = — (1 + n)/2 and its evolution in time is self- 
similar with logarithmic periodicity and with the same average law L(t) as for Burgers 
turbulence. Shocklines form self-similar regular tree-like structures. This model also 
describes important properties of the Burgers turbulence such as the self-preservation 
of the evolution of large scale structures in the presence of small scales perturbations. 
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1 Introduction 



The nonlinear diffusion equation 

dv dv d 2 v 
dt dx dx 2 



v(x, t = 0) = Vq(x). 



XI) 



was originally introduced by J. M. Burgers in || (1939) as a model for hydrodynamical tur- 
bulence. Burgers' equation (jl.lf ) describes two fundamental effects characteristic of any 
turbulence [IU]: the nonlinear redistribution of energy over the spectrum and the action 
of viscosity in small scales. Burgers' equation was used later to describe a large class of 
physical systems in which the nonlinearity is fairly weak (quadratic) and the dispersion is 
negligible compared to the linear damping |p9|| . The most important example of such waves 
are acoustical waves with finite amplitude ||24j| . Another class of problems, arising, e.g., in 
surface growth, also leads to Burgers' equation ||,||,[]30|. The three dimensional form of 
( |1 . 1|) has been used in cosmology to describe the formation of large scale structures of the 
Universe at a nonlinear stage of gravitational instability ( see e.g. Jl5|, [p5H , |T3| , |28|| ). 

In the physically important case of large Reynolds number, the action of viscosity is 
significant only in the small regions with high gradient of the velocity field. In the limit 
v — > 0, the solution of Burgers' equation has the following form (see 



my 



v(x, t) 



x - y(x,t) 



where y(x, t) is the coordinate of the maximum of the function 



G{x,y,t) = * (y) 



(X 



2t 



v {x) 



d^p(x) 
dx 



[1.2) 



1.3) 



Strong interaction between coherent harmonics leads to the appearance of local self- 
similar structures in Burgers' equation. A periodic initial perturbation with zero mean 
velocity is transformed asymptotically into a sawtooth wave with gradient d x v = 1/t and 
with the same period Iq. It is important that at this stage the amplitude a(t) = l^/t and 
the energy density o~ 2 {t) ~ /Q/12i 2 do not depend on the initial amplitude. 
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An initial one-signed pulse with the area m > 0, localised at t = in the neighbourhood 
of the point x = 0, also has asymptotically a universal form: it transforms into a triangular 
pulse with the gradient d x v — 1/t and increasing coordinate of the shock x s ~ (2mt) 1 / 2 . Due 
to the increase of the integral scale the amplitude of such a pulse a(t) = x s (t)/t ~ m 

l/2 t -l/2 

and its energy will decrease more slowly than for a periodic signal, like t^ 1 ^ 2 . 

Continuous random initial fields are also transformed into sequences of regions with the 
same gradient d x v = 1/t, but with random locations of the shocks separating them. Due 
to the multiple merging of the shocks the statistical properties of such random fields are 
also self-similar and may be characterised by the integral scale of the turbulence L{t). The 
merging of the shocks leads to an increase of the integral scale L(t), and because of this the 
energy 

a 2 (t) ~ L 2 {t)/t 2 (1.4) 

of a random wave decreases more slowly than the energy of periodic signals. 

The type of turbulence evolution is determined by the behaviour of the large scale part 
of the initial energy spectrum 

E (k) = a 2 k n b (k); (1.5) 
E (k) = ^J (v (x),v (x + z))e lkz dz. (1.6) 

Here bo(k) is a function which falls off rapidly for k > ko ~ Iq, and &o(0) = 1. For 
n > 1 the law of enrgy decay strongly depends on the statistical properties of the initial 



field (see e.g. [30] and references therein). For the initial Gaussian perturbation the integral 
scale Lit) ~ t 1 / 2 times logarithmic correction obtains and is determined by two integral 
characteristics of the initial spectrum: the variances of the initial potential ^/q an d the 



velocity v (x) 0,0, 0,0 

For n < 1 the structure function of the initial potential increases as a power law in space. 
Then the initial potential field is Brownian, or fractional Brownian motion, and some scaling 
may be used [f7|,Pm, [0' [EH> PLP^L p3| . In this case the turbulence is also self-similar 



and the integral scale L(t) increases as 

L(t) = (at) 2/(3+n) . (1.7) 
The energy of the turbulence is derived from (|1 ,4|) : 



At) ~ t -*, P= 2 ±±*L. (1.8) 
n + 3 

The difference between these two cases ( n < 1 and n > 1 ) is connected to the process of 
parametric generation of low frequency component of the spectrum. For the case n < 1 the 
newly generated low frequency components are relatively small and we have the conservation 
of large scale part of the spectrum: 

E(k,t) = E (k) = a 2 k n } for k «1/L(t). (1.9) 

Thus, the laws of turbulent decay are more complex than for simple signals, which can 
be attributed to multiple merging of the shocks. In QT2] a model of a regular fractal signal 
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with decay lower than for single one-signed pulse was introduced. The initial signal vq(x) 
was constructed as a sequence of one-signed pulses whose positions form a Cantor set with 
capacity (fractal dimension) D = In N/ In /3, where N p is the number of pulses in the scale 
L p ~ LiNf3 p ~ l , < D < 1. Multiple merging makes the decay of the wave slower and 
the general behaviour of the energy decay may be approximated by the power law with the 
exponent in (|1.8|) : 

P=\^, 0< P <1. (1.10) 

The evolution proves to be self-similar in successive time periods and (t i+ i,ti +2 ), 

where ti+i/tj = f3 2 /N. This shows log-periodical self-similarity of the field evolution. Linear 
and non-linear decay of fractal and spiral fields given by the sequences of regular pulses was 
also investigated in JT|. It was shown that the power law ( |1.8|) with the exponent given by 
the formula (|1.10|) holds also true for homogeneous fractal pulse signal with capacity D. 

Another model of a multiscale signal, which has the same general behaviour on the 
external scale L(t) (\L.7\j , and energy of the Burgers turbulence, was also discussed in [|12|]. It 



was assumed therein that the initial signal is a discrete set of modes - the spatial harmonics 



oc 



v o{ x ) = ^2a p sin(k p x + (p p ), (1.11) 

p=0 



with wavenumbers k p and amplitudes a p given in terms of a parameter e by 

k p = k e p , a p = a e~ hp , 



-(——), a = ak Q . (1.12) 



Amplitudes a p and the scaling exponent h are chosen from the condition that the mean 
energy of harmonics in the interval A p = k p — k p+ i be identical to that corresponding to the 
spectral density (|1.5| ): a 2 v = E(k p )A p . For e << 1 and n > — 1 the harmonics are spread 
over the spatial spectrum and accumulate at the point k = with decreasing amplitude. 
The main approach in this model was that the energy of the wave is the sum of energies of 
independent modes. The approach is nontrivial, but nevertheless leads to the same laws of 
the integral scale L(t) (|1.7|) and the energy decay ( |1.8|) as in the case of continuous spectrum 
(O). Let us point out that representation of the field given by the formula fll.llj ) is similar 



to shell models, which were introduced as useful models addressing the problem of analogous 
scaling in fully developed turbulence (see, e.g., |T1|, [13 and references in there). 



In present paper we consider the evolution of a regular signal whose behaviour in general 
is similar to the evolution of the Burgers turbulence with continuous spectrum (|1.5|) . The 
main difference, as compared to the model discussed above, is that we construct the exact 
solution of Burgers' equation using as an initial mode the "reverse" sawtooth wave. The 
frequency ratio in our model is e = l/N, where N is an integer and N > 2. These perturba- 
tion are similar to the well known Weierstrass and Weierstrass-Mandelbrot fractal functions 



sec 
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For the analysis of Burgers' equation it is convenient to use a mechanical interpretation. 
There is a one-to-one correspondence between the solution of Burgers' equation and the 
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dynamic of a gas of inelastically interacting particles [0, |14]]. Let us take a one-dimensional 



particle flux with a contact interaction: as long as the particles do not run into each other 
they move with constant velocity. In the collision they stick together, forming a delta- 
function singularity in the matter density. This leads to the appearance of gas of two 
species: a hydrodynamical flux of the "light" initial particles, and a gas of "heavy" particles 
arising in the adhesion process of light particles. The evolution of the particle velocity field 
will be described by the solution of Burgers' equation if we assume that the initial density 
of the light particles is po = const, the velocity of particles is equal to the initial velocity 
in fll. If) , and that the collision of the particles conserve their mass and momentum. This 
analogy permits construction of a very fast (linear time) algorithm of solution of Burgers' 
equation |27] . 



In our case, for the initial reverse sawtooth wave, all the matter turns into heavy particles 
at the same moment of time. Thus, after this time the evolution of the Burgers turbulence 
is fully determined by the motion of heavy particles, whose positions are positions of shocks, 
and masses are equal to Av ■ t, where Av are the amplitudes of the shocks. 

The paper is organized as follows. In Section 2 we consider the evolution and interac- 
tion of "reverse" sawtooth modes. In section 3 we consider the interaction of small scale 
mode with large scale structures. In Section 4 we investigate the properties of the sawtooth 
Weierstrass-Mandelbrot fractal function. In section 5 we show that deterministic model has 
logarithmic periodic self-similarity. We also discuss here the multi-dimensional generaliza- 
tion of this model. Section 5 presents concluding remarks. 

2 Evolution and interaction of "reverse" sawtooth modes 
in Burgers' equation 



Let us introduce the pth "reverse" saw tooth mode as 



V(p)(x, 0) = a p A(k p x + ip p ) 



(2-1) 



Here a p , k p are amplitude and wavenumber of the mode, — ip p is its phase. The function 
A(x) is 2tt periodic function, given on its first period by the following expression 



A{x) 



71 



x E [0,2tt[ 



(2.2) 



The set of "reverse sawtooth" functions is not orthogonal, but nevertheless we will introduce 
a set of modes satisfying equation (|1.12|) , whose wavenumbers and amplitudes satisfy the 

i.e. relation ( |1 . 1 1|) . We introduce the term 



same relations as the sinusoidal modes of [12 



"reverse sawtooth" because this signal is a sawtooth with teeth facing to the right, but the 
term "sawtooth" by itself is widely used in Burgers turbulence literature to refer to the late 
stage of the evolution of the wave profile a sequence of sawteeth with positive slope 1/t. 

The solution of Burgers' equation with a linear velocity profile Vq = —7(2 — x + ) is 
well-known (see, e.g., [p~3j| ) : 

-j(x-x+) 



v(x, t) 



1-7* 



(2.3) 
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The value 7 _1 has the dimension of time, and for 7 > 0, at the finite time t = 7 _1 the 
gradient d x v becames infinite. For 7 < the gradient becomes equal to d x v = t~ x , 
independent of the value of 7 at times t ^> |7| _1 . Thus, we have from ( |1.12j ), (p. 1[) , ( pT2|) , 
that the evolution of the pth mode is characterised by the nonlinear time 



t P = % l = a p k p ~ l = t /(e^y 2 r, t = lM n+3/2 (2-4) 

Based on solution ( |2.3| ) it is easy to see that for the "first" period (if <p p = 0) the evolution 
of the pth reverse mode at the initial stage (t < t p ) may be described as 

X A> ii0<x< f p V p ^ 

■ lp) (*,t) = { if lt _x| <lt ; (2-5) 



{x--f)lt, ifff<x-f <f . 

rZ<f) *-"rt n>w *^T) *^p 



On the other hand, at time t = t p , the mode transforms into a "direct" sawtooth wave with 
slope d x v = 1/t , independent of the amplitude and wavenumber of the mode: 

f x/t, if < x < jr; 

V(p)(x,t) = \r x _ if f < x < k . (2- 6 ) 

The density of energy a 2 {t) = (v 2 (x,t))L, where () L denotes averaging over the period, is 
conserved before t < t p , and decreases like {k p t)~ 2 after t > t p . 

Consider now the evolution of a gas of sticky particles in the case of independent evolution 
of the pth mode. In the general case, the density of the gas is calculated by using the Jacobian 
of the transformation from Lagrangian to Eulerian coordinates and may be written in the 
form ( see, e.g., ) 

p{x,t)=p Q {l-td x v{x,t)) (2.7) 
Then it is obvious that at the initial stage 

. . 1 . tx t ix . . 

p(x,t)=po , |_-a;|<-— (2.8) 

while p is zero outside this interval in each period. At time t = t p all the light particles in 
each period collide into a single heavy particle with mass 

2tt . . 

m p = p L p = p — , (2.9) 

K p 

and the heavy particles have positions 

* p ,z = f-^ + ^Z; Z = 0,±l,±2,... (2.10) 

equal to the zero positions of the initial pth mode. The process of light merging particles 
and the evolution of the velocity is shown in Fig.l. 
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Consider now the joint evolution of two successive modes: pth and (p+ l)th. From (|2.4|) 
one can see that the ratio of nonlinear times of the successive modes is 



W// ~~ tp,p+i — 7 — ~. — ~~ 1 — , I P u ' (2-12) 



W = e -(n+3)/2 = e l-fc (2 n) 

which does not depend on p and increases if the exponent n is greater than —3. The gradient 
of the initial field %>)(#) + f( p+ i)(x) is — (j p + 7 P +i), so the effective nonlinear time for 
such a sum is 

1 = t ± 

7p + 7p+i 1 + t p / tp+i ' 

Because all parts of the initial perturbation have the same slope, all light particles will collide 
at the same time t = t PjP+ i. 

The mass of heavy particles after merging are m Pti = p Aj )i+1 , where A^+i is the distance 
between adjacent shocks in the initial perturbation. The ratio of periods of two adjacent 
modes is L p+1 /L p = k p /k p+ i = e" 1 . UN — e~ x is an integer larger than 1, there will be 
N + 1 heavy particles on the period of the larger scale mode (p + l)th: (N — 1) with the 
mass m p = poL p (|2.9|) , and two particles with total mass equal to m p . These two particles 
only exist when the shock of the (p + l)th mode is located in the interval between shocks of 
the pth mode. For simplicity, we will consider the case where the spatial relations 

k p+1 ip p+1 = k p ip p + 2-nr/N, (2.13) 

between the phases of successive modes, hold. In this case the discontinuities of the (p+ l)th 
mode do not produce new shocks in the total perturbation Vi^(x)+Vuj + i)(x). Thus, at times 
t larger than t > t r>e ff, the masses of all heavy particles will be the same, as would be the 
case without the large scale modes (|2.9| ). 

The positions of these heavy particles at time t = t Pt£ ff are 

X(p,l)(tp, e ff) = Xp,l + V(p + i)(x p! t)t Pte ff, (2-14) 



where x P) i are the zero positions of the pth. mode ( |2.1(J| ). The velocity of this particle is equal 



to V(p + i)(x Pt i). Equation fl2.14|) is obvious if we use the trivial equality V( p )(x P) i)+V( p+ i)(x P} i) 



V( p+ i)(x p j), and also note that the position of the heavy particle x P) i{t P)e jf) is equal to the 
position at the same time of all light particles with initial coordinate x = x p j. From ( |2.14| ) 
we immediately have that after time £ p , e // the positions of the particles are 

X(p ; i) (t) = x {p j) + V(p + i) {x P; i)t. (2.15) 

The difference between the coordinates of the adjacent particles X P) i{t) and X p j +1 (t) de- 
creases with time, proportionally to the gradient of f( p+ i)(x): 

dv (^*) 

(t) — X Pi i(t) = (x p j + i — x Pi i) — t — {x P ,i+i — % P ,i) = 

(x Pj i +1 -x Pj i)(l-t/t p+1 ). (2.16) 
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These particles collide at time t = t p+ \ (|2.4j) and the newly created heavy particles will 
have masses 

m {p+1) = p L p+1 (2.17) 



and positions 



x p+ljl =- + - ; / = 0,±1,±2,... (2.18) 



The velocity of these particles is zero. 

Thus, at times t larger then t p+ i, the evolution of the initial perturbation vo(x) = v ( p )(x) + 
v (p+i) ( x ) will be the same as the evolution of only the large scale mode (x). The process 
of particles merging and the evolution of the velocity for the sum of two successive modes 
with the periods ratio N = 2 are shown in Fig. 2. 

By recurrence, it is evident that for finite number of modes v(x) = %>)(#) + t>( p+ i)(x) + 
... + V(m)(x) the evolution of the field after tu will be the same as the evolution of only 
the largest mode V(m\{x). The reason for this, is of course, the special relation between 
the phases <p p and the wavenumbers k p of all interacting modes: k p = ko/N p ( see equation 
( |1.12| ) ). For integer N the minimal value of any combination of these wavenumbers is 
equal to the largest mode wavenumber ku = ko/N AI . So the nonlinear interaction does not 
produce new components at frequencies less than ku- 

3 Interaction of small scale "reverse" sawtooth mode 
with large scale structures 



Let us now consider the interaction of the pth mode with an infinite series of larger scale 
modes 

OO CO 

W p (x) — v r (x) = a r A[k r x + (p r ), (3.1) 

r=p+l r=p+l 

assuming that the phases of the modes satisfy the relations ( |2.13j ) and that k r = koe r , a r = 
a e~ hr , h = — (n + l)/2. From ( [3.1|) and (2^) we have for the gradient of the initial 
perturbation v (x) = v p (x) + W p (x) 

d x v (x) = d x v p (x) + d x W p {x) = J2™= p lp = 

= 7oE r =p(e 2 ) =7oe 2 p — prrsr, (3.2) 

the condition n > —3 {h < 1) being necessary for the series to converge. From ( |3.2|) , we 
have for the effective time of nonlinearity of pth mode 

(n+3) 

t p = l/d x v {x) = t p {l-e—), (3.3) 



with the original t p determined by the equation (p.4j ). 
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Thus, after the time of collision t p , heavy particles with mass m p = poL p ( |3.2|) appear. 
The coordinates of these particles will be determined by an equation similar to (|2.15 ) 



x 



P M = + W p (x Pil )t, (3.4) 



with the velocity of particles determined by the function W p (x) (see |3.1| ), which is a sum of 
all larger modes, and the coordinates of the zeros of the pth modes. The difference 

between the coordinates of adjacent particles x p j and x p j will then decrease with time like 
(1 — t/tp+i), where = £ p +i(l — e^ n+3 ^ 2 ) is the inverse of the gradient of the function 
W p (x), see (|3.1|) and (3.2). Thus, the time of particle collision for this generation will be 



described by equation (|3.3| ) with p — p + 1, and the new masses will be determined by the 
period of the (p + l)th mode - see equation (|2.17|) . 



The extrapolation of this particle merging process to the next generations is evident by 

.— — ~ ti-\-3 

recurrence. The qth collision of heavy particles takes place at time t p+q = t p+q (l — e~*~), 
the masses of these particles at this time are determined by the period of the (p + g)th mode 
m p+q = p L p+q = 27cpo/k p+q (P0|), ( |2.17p . In the time interval t G [t p+q ,t p+q+ i], 



tp+q+l (n+3) , 

t p +q 



2 



N^- (3.5) 



the coordinates of particles will be determined by the equation ( |3.4| ) with p = p + q. Here 
W p+q (x) is the sum of the velocities of all larger modes with r > p + q, and x p+qt i are the 
zeros of (p + g)-th mode. It is important to note that at time t > t p+q the evolution of the 
particles is solely determined by the modes with r > p + q. It means, that at times t > i p+q 
the position of the particles does not depend on the presence in the initial condition of the 
small scale modes with r < p + q. 

Thus, two processes with different initial velocities: Vo(x), the field with small scales, 
and vo(x), the field without small scales modes: 

v (x) = W p+q -x(x)] v (x) = W p -i(x) (3.6) 

will have the same evolution after t > t p+q -\. Even if p — ► — oo (when modes with very 
small scales L p ~ e _p = N p and very large amplitudes a p ~ ao(e^ n+1 ^ 2 ) p = a (N ( - n+1 ^ 2 )^ p ) 
are present in the initial perturbation) the multiple merging of the particles will lead to the 
independence of the evolution of large scale modes with respect to the small scale modes. 
This effect is similar to the self-preservation of large scale structures in Burgers turbu- 



lence ||,[|T7|]. When the initial field vo(x) is noise, the highly nonlinear structures contin- 
uously interact and due to the merging of shocks, their characteristic scale L(t) constantly 
increases. The presence of small scale noise perturbation Vh{x) results in additional fluc- 
tuations in the shock coordinates Axk(t), and these fluctuations increase in strength with 
the passage of time. Thus, the final result of the evolution of the field is determined by the 
competition of two factors, the increase in the external scale L(t) of the structures and the 
increase in the strength Axk(t) of shock coordinates fluctuations, the later being related to 



the perturbation Vh(x). In a turbulence, having power index n < 1 (1^5), multiple merging 



of shocks leads to self-preservation of the large scale structures independently of the presence 
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of small scale components. For the model signal this effect appears for arbitrary n due to 
the special choice of wavenumbers and phases of interacting modes. 

It was stressed in the introduction that the solution of Burgers' equation has a one-to-one 



correspondence with the dynamics of the gas of inelastically interacting particles ( p3|1 ). The 
stage when all light particles collide, forming heavy particles, corresponds to the solution 
of Burgers' equation with a well-defined slope d x v = 1/t. In this case the profile of the 
field v(x,t) is fully determined by the coordinates and amplitudes of the shocks. Their 
coordinates X s (t) are equal to the coordinates of heavy particles, their velocity 

v s (t) = = {v s {x s -0,t) + v s {x s + 0, t))/2 (3.7) 

is equal to the velocity of the particles, and the amplitude of the shock 

Av s (x) = (v(x s - 0, t) - v(x s + 0, t)) = m/t (3.8) 

is determined by the mass of the particle (po = 1) (see, e.g., |H§). 

Thus, the investigation of the motion of heavy particles permits to fully reconstruct the 
properties of the velocity field v(x, t) of Burgers' equation. 



The sawtooth Weierstrass-Mandelbrot fractal func- 
tion 



It was shown in the previous section that the evolution of the particles (shocks) is deter- 
mined by the function W p (x) (|3.1| ). The basis functions of W p (x) are the reverse sawtooth 
periodic functions with wavenumbers k r = k e r and amplitudes CL r — Ojq(. ~ h \h= -(n + l)/2 
satisfy relations ( |1.12| ). Wavenumbers form a geometrical progression like in the Weier- 
strass function (see [|J) and accumulate at the origin k = 0. In the original Weierstrass 
function, the situation was the opposite with increasing frequencies, but nevertheless the 
function W p (x) has many properties of Weierstrass function and of its generalisation - the 
Weierstrass-Mandelbrot function (see pl| , [|J). 



We consider here a deterministic function W p (x) with the special phase relation <p p = 
(2itk/N)p (k = 1,2, N; N = 1/e ), thus, the discontinuities in the largest modes 
r > p + 1 coincide with some of the discontinuities of the smaller mode r = p + 1. The 
function W p (x) is continuous in the intervals 27r/fc p+ i = 2n/{koe p+1 ) with the same slope in 
each interval. The inverse value of this slope t p+ i 

£ p+1 = t p+1 (l-e^); t p+1 = t (e^) p+1 (4.1) 

is proportional to the nonlinear time t p+ i of the smallest mode. Of course, we need n > — 3, 
so that the convergence of ( ^.2[ ) is assured and the inequalities t p+ i > t p hold. The amplitudes 
of the modes are proportional to e^ n+1 ^ 2 and for n > — 1 the function W p (x) is bounded 

OO 1 

W p (x) < £ a r = a (e^r +1 ^ . (4.2) 

r=p+l 1 — e 2 
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Thus, for finite p the energy of W p (x) is also finite. For the case of the phase relation 
introduced above, the functions W p (x) also have scaling properties, so that for instance for 
k = 0, we have 

W p (x) = e- hp W (e p x); W p (e m x) = e hm W m+p (x). (4.3) 

The case — 1 < n < 1 is similar to the initial conditions with generalized white noise 
in Burgers turbulence. The energy of the initial signal in such turbulence is determined by 
the largest cutoff wavenumber, so in our model by the smallest scale p. If p — > — oo the 
energy of the model signal (as the energy of white noise) will tend to infinity. But from 
the considerations in the previous section we have that at the finite time t all the modes 
with t p < t have finite energy ~ L 2 /t 2 due to the nonlinear dissipation, so that the whole 
energy of the turbulence is also finite. Thus, even in the case of " divergent" initial conditions 
(p — > — oo), we will have a "convergent" solution for any time t > 0. 

The case n < — 1 is similar to having fractional Brownian motion initial condition in 
Burgers turbulence. In this case, the series ( [4.2[ ) diverges and the initial signal W p (x) is 
unbounded. But for Burgers turbulence ( for the process of particles motion and colli- 
sions ) only relative velocity of the particles matters. So we can use the same regularisation 
procedure with W p (x). Such a procedure was done with the Weierstrass function in 

In our case, taking <p p = 0, we can introduce the function W£°(x) = W p (x) — W p (0), 
according to which is finite in all finite spatial intervals. The other way to get a bounded 
function is to use special phase relations for the modes. 



5 Self-similarity properties of deterministic model in 
one and two dimension 

Here we summarise the properties of the evolution of the multiscale deterministic signal using 
some additional information about scaling characteristics of W p (x), and compare them with 
the properties of the Burgers turbulence. 

Let us consider the evolution of the multiscale signal 

v Q (x) = v p (x) + W p (x). (5.1) 

It was shown that at times t for which t > t p , heavy particles with mass M p = poL p appear 



and their coordinates are determined by the relation (|3.5| ). These particles collide at time 



t p +i, (t p+ i/t p = N ( 2 5 ; TV = e x ), and new particles with masses m p+1 = p L p+1 = m p N 



appear. Their motion will be determined by the same law ( |3.5| ) with substitution p — > p + 1. 
Using the scaling properties of W p (x) fl4.3| ) we have, that the motion of the particles in this 



interval will be similar to the motion of the particles in the interval t p+ i] if we rescale 
the time t/t p =>- t/t p+ \. Since the ratio t p+ i/t p does not depend on p, one can speak about 
the logarithmic periodic self-similarity of the motion of the particles. This means that at 
arbitrary interval [t q , t q+ i] the motion of the particles will be similar to the motion of the 
particles in the interval i p ,t p+ i, by the scaling factor x p /x q = e p ~ q in space, and the scaling 

n _i_3 

function t p /t q = [e~^~) p ~ q in time. The coordinates and masses of the particles fully 
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determine the velocity field, and so the solution of Burgers' equation is also logarithmic 
periodic self-similar. 

With each collision, the mass M(t) of the particles increases N = 1/e times. The time 

n _i_3 

interval between the two successive collisions increases as t p+ i/t p = N^~. Thus, by the 
approximation of piecewise constant function m(t) by the power law 

m(t) ~ m (t/t ) (n+3)/2 - (5.2) 

we obtain the same result as for the Burgers turbulence. In our case, m is proportional to 
the period of the smallest mode at time t, and is analogous to the integral scale in Burgers 
turbulence. 

In the case n > — 1 we can also estimate the energy decay of the model signal. For 
n > — 1 and e < 1 the main energy of the signal at time t is in the smallest mode and is 
proportional to L 2 (t)/t 2 . Thus, we have here again the same law for the energy decay as for 
Burgers turbulence. 

The numerical simulation based on the algorithm [[27|] was done to illustrate the process of 



particles merging and velocity field evolution. The trajectories of the particles and profile of 
the field at different times are plotted for the initial "white noise" signal (n = 0, h = —1/2) 
in Fig. 3, and for the initial "Brownian" motion (n = —2, h = 1/2) in Fig. 4. Ten modes with 
the ratio of successive wavenumbers e = 1/N = 1/2 were used. The plots show the initial 
stage of the evolution in some relatively small region where the finiteness of the number of 
modes is not significant. 

In Fig. 3 one can see that for n = the initial "sawtooth" multiscale function oscillates 
near v — like a "white noise" with finite variance. After the collision of light particles, 
when the reverse sawtooth function transforms into a sawtooth wave with positive gradient 
d x v = l/t, the structure of the signal is relatively simple, and even for N = 2 the main 
energy remains in the mode with smallest wavenumber. 

In the case n = —2 the initial profile has a large deviation behaviour which is typical for 
Brownian motion functions. After merging of light particles, the sawtooth profile has a set 
of small shocks with different amplitudes, which is also similar to the properties of Brownian 



signal in the Burgers turbulence 28 



In Figs. 3(c) and 4(c) the velocity field at three successive merger times t*/t** = N^ n+3 ^ 2 
are plotted. These figures show the logarithmic periodic self-similarity of the evolution of 
multiscale signals. 

We notice now that such multiscale waves may be constructed for multidimensional 
Burgers' equation. Let us assume that the initial vector field V p (x) is an infinite series of 
"reverse" modes v r (x): 

oo 

V p (x) = £ v r (x) , (5.3) 

r—p 

In the two dimensional case, the rth "reverse" mode may be composed of piecewise linear 
functions defined on a system of regular triangles of size L r covering the plane. We consider 
here the special case when the ratio of the scales of two adjacent modes is L r+ i/L r = e _1 = 
N = 2. We assume also the special symmetry and phase relation between the different 
modes. In our case one big triangle is divided into four smaller triangles with vertices 
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located at midpoints of its sides ( see Fig. 5 ). We assume that inside each triangle in the rth 
mode the velocity has a linear profile v r (x) = — 7 r (x — x + ), where x + is the coordinate of 
the center of the triangle. The solution of the multidimensional Burgers equation for such 
initial perturbation 

v(x, t )= -^ X 7+> (5.4) 

1 — jt 

is now valid inside the triangle of size L r (t) = L r (l — tj r ) . The value 7" 1 has the dimension 
of time, and at the finite time t r = 7" 1 the velocity gradient becomes infinite. 

On the other hand, at time t = t r , the mode transforms into a "direct" sawtooth wave 
with the universal behaviour inside the new set of triangles and with the gradient 1/t inde- 
pendent of the amplitude and wavenumber of the mode: 

v(x,i) = ^p, (5.5) 

where x c is now the center of the triangle, coinciding with the top of the initial triangular set. 
Consider now the evolution of a gas of sticky particles in the case of independent evolution 
of the rth mode. Then, it is obvious that at the initial stage, inside the "collapsing" triangle 
of size L r (t) = L r {l — tj r ) the density increases as 

p(x, t) = po-. 777^2 > (5-6) 

(l - t/t p y 

while p is zero outside "collapsing" triangular in each cell. At time t = t r , all the light 
particles in each cell collide into a single heavy particle with mass 



p ^v / 3/4, (5.7) 



and the heavy particles' positions are equal to the center of the initial triangle x + . 

We assume also that the evolution of the rth mode is characterised by a nonlinear time 
t r the same as in the one dimension case 



U = i; 1 = V(2- (n+3)/2 ) r . (5.8) 



Let us now consider the evolution of the vector field V p (x) ( |5.3|) which is an infinite 
series of "reverse" modes. The evolution of the vector field is very similar to the evolution 
of the scalar field (|3.1|) . For the gradient of the initial perturbation V p (x) we have the same 
relation ( |3.2| ) as in ID. The effective time of nonlinearity of the smallest pth mode in the 
vector field (|5.3| ), in presence of all large scale modes, is determined by the equation (|3.3|) . 
Thus, after the time of collision t p , heavy particles with mass m p Q5.7D appear. Velocities of 
these particles will be determined by the function V p+ i (x) ( |5.3| ), which is a sum of all larger 
modes, but the number of particle collisions is determined by the next p + 1 mode. At time 
t p+ i we have a collision of four heavy particles. 

The extrapolation of this particle merger process to next generations is evident by recur- 
rence. The gth collision of heavy particles takes place at time t p+q ), ( |5.8| ) the masses of these 
particles at this time are determined by the scale of the (p + g)th mode m p+q = p Lp +(? a/3/4 
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(|5.7| ). Here also one can speak about the logarithmic periodic self-similarity of the motion 
of the particles. This means that at arbitrary interval [t q , t q+ i\ the motion of the particles 
will be similar to the motion of the particles in the interval t p ,t p+ i, by the scaling factor 
Xp/x q = 2^~ p ) in space, and the scaling function t p /t q = (2^)^^ in time. 

We used computer simulation for studying two dimensional case; results of the simulation 
were generated in so called VRML (Virtual Reality Modeling Language), which enables to 
handle three dimensional figure in different projections. Fig. 6 presents snapshots of this 
modeling. On the Fig. 6(a) one can see that particles formed by small triangles move towards 
the center of an embracing triangle; this center in its turn, moves towards the center of the 
next bigger triangle in the hierarchy; e.t.c. Fig.6(b) gives the side-view of this process. 

6 Conclusion 

In conclusion, we would like to point out that the evolution of the multiscale signal with 
the Weierstrass spectrum simulates properties of Burgers turbulence such as self-similarity, 
conservation of large scale structures and has the same laws of the energy decay and integral 
scale. The difference between the deterministic model and Burgers turbulence is that here 
we have the exact solution for the evolution of multiscale signals and these properties are not 
stochastic but deterministic. The evolution of the multiscale signal is exactly self-similar in 
logarithmically spaced time intervals. The evolution of the large scale modes is completely 
independent of the small scales modes, even if these have very large amplitudes. 

These properties take place for Burgers turbulence in the stochastical sense and, more- 
over, for a signal with cutoff frequencies of small scales, only asymptotically. Of course, these 
properties of the multiscale signal are determined by the special form of modes (reverse saw- 
tooth function), the special relations between wavenumbers of modes (Av+i = k r /N, where 
iV is an integer) and their phase relations. 

On the other hand, these model signals do not reflect such properties of Burgers turbu- 
lence as qualitative difference in the behaviour of the turbulence for n < 1 and n > 1 in 
the power spectrum ( |1.9| ) due to the process of generation of large scale components in the 
spectrum. For the deterministic model this process is not present due to the special relation 
of wavenumbers. 

Let as now move to the mechanical interpretation of solution of the Burgers equation. 
For the initial reverse sawtooth wave, all the matter turns into heavy particles at the same 
moment of time. Thus, after this time, the evolution of Burgers turbulence is fully deter- 
mined by the motion of heavy particles. The trajectories of heavy particles form regular 
tree-like structure on the plane (X, t), see Figs. 3 and 4. The properties of this structure 
depend on the parameters of our model. The integer e _1 = N is the number of trajectories 
which intersect at one point and form a new branch of our structure. For e _1 = N = 2 we, 
thus, obtain binary tree structure. Changing of the parameter h stretches or contracts the 
structure in the t direction. One can say that our structure is the plane representation of 
the N— tree; the root of our tree is located at t — +oo. This tree is similar to the flattened 
fractal model of botanical umbrella tree (see pifl ). If we take some node (X, T) of this 
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structure as a root, and consider the trajectories of all heavy particles, which will merge 
at the moment of time T at this point (X, T) we shall also obtain an n-tree. The whole 
tree seems self-similar, because every branch plus the branches it carries is a reduced scale 
version of the whole. 

We also notice that such multiscale waves may be constructed for multidimensional 
Burgers' equation. 
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Figure captions 



Figure 1: Evolution of one-mode pulse, (a) particle trajectories; (b) evolution of the initial 
velocity field given in the same spatio-temporal scale. Bold lines on the time axis denote 
moments at which the profiles of the velocity are plotted. 

Figure 2: Evolution and interaction of two modes, (a) particle trajectories; (b) evolution of 
the initial velocity field given in the same spatio-temporal scale. Bold points on the time 
axis denote moments at which the profiles of the velocity are plotted. 

Figure 3: Evolution of the multiscale fractal signal with n = (h = — 1/2), corresponding to 
"white noise" signal, (a) particle trajectories; (b) evolution of the initial velocity field; (b) 
velocity field taken at the initial moment of time and then at three successive time moments 
of self-similarity. 

Figure 4: Evolution of the multiscale fractal signal with n = —2 (h = 1/2), corresponding 
to "Brownian motion" signal, (a) particle trajectories; (b) evolution of the initial velocity 
field; (b) velocity field taken at the initial moment of time and then at three successive time 
moments of self-similarity. 

Figure 5: Plane construction for two dimensional case. The hierarchy of triangles, used 
for the construction of the multiscale signal, with four layers shown. The initial signal is 
constructed as a series of signals piece-wise linear on triangles. 

Figure 6: Particle trajectories for the multiscale fractal signal in two dimensional case pre- 
sented in spatio-time three dimensional space; the width of particle trajectory reflects its 
mass: (a) top view; (b) side-view. 
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